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ABSTRACT 


The  scale  height  in  height-bands  between  150  and  400  km  is  assumed 
to  vary  linearly  with  height.  Integration  of  the  hydrostatic  equation  for  an 
ideal  gas  above  a  spherical  earth  then  leads  to  a  power  function  represent¬ 
ation  of  the  air  density  over  the  band.  With  integral  exponents  such  power 
laws  give  better  fits  to  several  proposed  model  atmospheres  over  altitude 
ranges  of  several  hundred  kilometers  than  those  provided  by  the  usual 
exponential  representation  of  air  density. 

The  representation  of  air  density  in  satellite  deceleration  equations 
by  power  functions  with  integral  exponents  reduces  them  to  elementary 
forms.  It  has  been  possible  with  such  density  distributions  to  obtain 
simplified  formulas  which  may  be  useful  for  (a)  computing  atmospheric 
densities  from  satellite  accelerations,  (b)  comparing  proposed  model 
atmospheres  with  observations,  and  (c)  developing  further  the  theory  of 
satellite  orbits  in  the  presence  of  air  drag.  As  is  possible  with  the  expo¬ 
nential  form,  these  power  functions  may  be  modified  to  take  account  of 
the  effect  of  an  oblate,  rotating  atmosphere.  Their  use  may,  therefore, 
permit  the  development  of  a  simplified,  accurate  orbit  theory  for  satel¬ 
lites  with  perigee  heights  below  300  km.  Certain  preliminary  results 
are  discussed  and  compared  with  previous  theory  and  observations. 
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INTRODUCTION 


The  orbit  of  a  close  earth  satellite  is  continuously  contracting. 

This  is  because  the  main  effects  of  air  drag  are  to  reduce  the  eccentri¬ 
city  e  and  the  length  of  the  semimajor  axis  a.  The  contraction,  which 
is  most  pronounced  when  the  perigee  height  is  below  300  km,  is  well 
described  by  two  equations  [King-Hele,  Cook  and  Walker,  1959;  Parkyn, 
i960].  These  equations  express  the  changes  Aa  and  Ae  in  the  semimajor 
axis  and  eccentricity  during  one  revolution  from  perigee  to  perigee  in 
the  form 


2  Ztt 

Aa  =  -  6a  f  (1 

••o 

2 


Ae  =  -  6a  (1  -  e  ) 


3/2  -1/2 

+  e  cos  E)  {1  "  e  cos  E)  pdE 


(1) 


27r  1/2 

cos  E  (1  +  e  cos  E)  (1 


e  cos 


pdE 

(2) 


where 


6  =  FCj^  A/M 

F  =  factor  to  account  for  rotation  of  atmosphere 


Cjj  =  drag  coefficient 
A  =  effective  cross  section  of  satellite 


m  =  mass  of  satellite 


E  =  eccentric  anomaly 
p  =  atmospheric  density. 

Equation  (1)  has  been  employed  in  the  past  as  a  basis  for  the  devel¬ 
opment  of  formulas  for  computing  atmospheric  densities  from  observa¬ 
tions  of  satellite  decelerations  [Sterne,  1958;  Groves,  1958;  King-Hele, 
Cook  and  Walker,  1959].  When  coupled  with  (2),  (1)  also  provides  a  log¬ 
ical  starting  point  for  the  development  of  equations  which  relate  the  vari¬ 
ations  of  a  and  e  with  each  other,  with  the  time,  and  with  the  revolution 
number.  To  make  use  of  (1)  and  (2),  an  explicit  functional  dependence 
of  air  density  on  the  height  r  -  q  above  a  reference  level  q  must  be 
assumed.  The  reference  level  and  r  are  to  be  measured from  the 
geocenter.  ® 
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In  many  previous  investigations  [e.  g.  -  King-Hele,  Cook  and  Walker, 
1959:  Parkyn,  I960]  an  exponential  form  of  the  vertical  distribution  of  air 
density  ,(r)  e  p(,^)  enp  (  -  ,r  -  has  been  assumed  in  »hich  Ute 

density  scale  height  H'  is  treated  as  constant  in  the  integrals  of  (1)  and 
(2).  Evaluation  of  these  satellite  deceleration  integrals  then  leads  to 
formulas  which  involve  modified  Bessel  functions  of  the  first  kind  of 
argument  ae/H'.  If,  as  Jacchia  has  pointed  out  [Jacchia,  I960],  H'  ex¬ 
hibits  an  altitude  dependence,  a  systematic  error  is  introduced  into  the 
development. 

An  examination  of  the  analyses  and  data  presented  in  Nicolet[Nicolet, 
1959]  and  several  preliminary  or  proposed  model  atmospheres  [Kallmann, 
1959;  Minzner,  Champion,  and  Pond,  1959:  Paetzold,  1961]  suggests  a 
linear  fit  to  a  given  model  as  a  first  approximation  to  the  variation  of  the 
atmospheric  scale  height  H  with  altitude  in  a  height  band  of  interest  within 
the  altitude  range  considered  here  150-400  km.  These  height  bands  can 
be  fairly  broad,  covering  a  range  of  from  100  to  200  km.  For  example, 
in  the  height  band  200-400 'km,  of  greatest  importance  in  satellite  orbit 
theory  in  the  presence  of  air  drag,  the  scale  height  over  the  entire  band 
may,  for  the  Nicolet  data  and  for  two  of  the  models  cited,  be  reasonably 
well  approximated  by  constant  gradients.  The  third  model,  due  to  Paetzold, 
can  be  approximated  by  a  linear  fit  over  the  height  band  220-400  km. 

In  the  discussion  in  the  following  section  it  is  seen  that  the  assump¬ 
tion  of  a  linear  variation  of  scale  height  H  over  a  given  height  band  leads, 
through  the  equation  of  hydrostatic  equilibrium,  to  a  modified  power  func¬ 
tion  representation  of  air  density  over  the  band.  Similarly,  the  assump¬ 
tion  of  a  linear  variation  of  H'  with  altitude  leads  to  a  pure  power  function 
representation.  In  either  form,  the  exponent  n  of  the  power  law  depends 
only  upon  the  value  of  the  constant  scale  height  gradient  and,  in  general, 
it  will  not  turn  out  to  be  an  integer.  The  introduction  of  such  a  vertical 
distribution  of  atmospheric  density  into  (1)  and  (2)  leads  to  integrands  for 
which  no  simple  integration  is  apparent.  On  the  other  hand,  power  func¬ 
tions  with  integral  exponents  reduce  the  satellite  deceleration  integrals 
to  elementary  forms  and,  hence,  may  provide  a  base  for  developing  form¬ 
ulas  which  may  be  useful  for  (a)  computing  atmospheric  densities  from 
satellite  accelerations,  (b)  comparing  proposed  model  atmospheres  with 
observations,  and  (c)  developing  further  a  simplified  satellite  orbit  theory 
in  the  presence  of  air  drag.  The  possibility  of  so  representing  the  air 
density  over  the  interval  150-400  km  will  be  considered,  and  certain  pre¬ 
liminary  results  derived  from  the  use  of  such  a  representation  will  be 
introduced  and  discussed. 
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ANALYTICAL  REPRESENTATIONS  OF  AIR  DENSITY 
BASED  ON  A  VARIABLE  SCALE  HEIGHT 


Values  of  atmospheric  density  and  scale  height  taken  from  Nicolet, 
[Nicolet,  1959]  and  the  several  model  atmospheres  previously  cited  are  pre¬ 
sented  in  Table  1  for  the  altitude  interval  1 50-400  km.  The  same  data  are 
showninFigs,  landZ.  In  Fig.  1  itis  seen  that,  except  for  the  band  1 80-220 
km  in  the  Paetzold  data,  there  is  a  significant  steepening  of  the  slope  d  log  p 
/dr  below  300  km  which  increases  (negatively)  as  the  altitude  decreases. 
Thus,  an  exponential  distribution  of  air  density,  which  would  appear  as 
a  straight  line  on  Fig.  1,  could  not  reasonably  be  expected  to  represent 
the  data  well  below  300  km  except  over  limited  intervals. 

The  curves  depicted  in  Fig.  2  have  been  constructed  by  the  rather 
arbitrary  procedure  of  connecting  the  published  data  points  with  straight 
lines.  Thus,  they  indicate  only  the  gross  character  of  the  data.  Reason¬ 
able  linear  approximations  to  these  data  may  be  obtained  by  treating  the 
two  height  bands  150-200  km  and  200-400  km  separately  for  all  but  the 
Paetzold  model.  For  each  of  the  two  bands  and  for  each  set  of  data,  a 
gradient  of  scale  height  may  be  assigned.  For  the  Paetzold  data  three 
bands  150-180,  180-220,  and  220-400  km  may  be  treated  similarly. 


Assume  then  a  linear  variation  in  scale  height 


H  =  H  +  /3  (r  -  q  ) 
0  0 


(3) 


where  p  is  the  scale  height  gradient  and  Hq  is  the  value  of  H  at  the  refer¬ 
ence  level  q  .  Let  the  intercept 


s  =  q  -  H  /  (4) 

0  0 

of  the  line  (3)  with  the  r-axis  be  introduced  into  (3)  and  the  result  substi¬ 
tuted  into  the  equation  of  hydrostatic  equilibrium  adapted  to  a  terrestrial 
gas  [Nicolet,  1954] 


Hdp  =  -  pdr. 


(5) 


Equation  (5)  may  be  integrated  to  yield 

P  =  Pq  [(qg  -  s)/(r  -  s)]  (6) 

where  p^  is  the  value  of  the  atmospheric  pressure  p  at  the  reference 
level  q^  and  n  =  1  +  P  .  The  air  density  may  be  introduced  into  (6) 

from  the  perfect  gas  law 
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Pi 

w 

^  * 
50 


1 


5 


Fig.  1.  Variation  of  p  with  Height  Above  the  Earth's  Surface  for  Various  Atmospheric  Models. 


where  is  the  value  of  the  acceleration  of  gravity  g  at  the  reference 

level  q  .  This  substitution  yields  the  modified  power  function  for  the 

0 

air  density 

(2)  2  n 

P  =  Pq  -  s)  ]  .  (7) 

The  superscript  notation  on  the  symbol  for  air  density  is  introduced  into 
(7),  and  elsewhere,  for  convenience  in  distinguishing  among  the  different 
air  density  functions  to  be  discussed.  For  the  special  case  that  H  is 
constant  (7)  reduces  to  the  modified  exponential  form 

(2)  2 

lim  p  = 

iS  -0 

Alternative  forms  of  (7)  suggest  themselves  since  the  factor 

is  close  to  unity  in  the  region  of  interest.  For  example,  the  functions 

(1)  n 

P  =  pQ  -  ®)  /(J*  -  s)l 


(0)  n 

P  =  P^[  (q^  -  s)/(r  -s)]  (9) 

give  values  of  air  density  which  differ,  in  general,  by  less  than  three 
per  cent  over  the  entire  band  200-400  km  from  those  ,'i'iven  by  (7). 


The  form  of  (8)  will  find  application  later  in  developing  a  relation¬ 
ship  between  the  semimajor  axis  and  the  eccentricity.  The  form  of  (9) 
may  be  derived  from  the  definition  of  density  scale  height  [Jacchia,  1959; 
King-Hele,  Cook  and  Walker,  I960] 

H'dp  =  -pdr.  (10) 


If  a  linear  variation  of  density  scale  height  with  altitude  H'=  /3'(r  -  s') 

is  assumed,  integration  of  (10)  yields  p^^^  =  p^  [(q^  -  s')/(r  -  s')]^  . 

Here  s'  =  q^  -  /  /3  ',  is  the  density  scale  height  at  q^  and  ft  =  1/n' 

is  the  density  scale  height  gradient.  If,  however,  H'  is  taken  to  be  con¬ 
stant  (10)  leads  to  the  more  familiar  exponential  representation  of  air 
density 


(3) 

P 


(0) 

=  lim  p  =  p^  exp  [  -(r  -  qp)/H']. 

ft  '-0 


(11) 
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If  a  linear  variation  of  y  ■  Hg  is  assumed  another  air  density  dis¬ 
tribution  which  combines  features  of  the  power  and  exponential  laws  may 
be  obtained.  Take 


y  =  Vq  +  K  (r  -  q^)  =  K  (r  -  S) 


(12) 


where  Y  is  the  value  of  y  at  q  ,  K  is  the  Y  gradient,  and  S  =  q«  -  Y ^IK. 

0  0  u  u 

Substitution  of  (12)  into  (5)  leads  to 

^  =  Pq  (r/qg)’^  [  (qQ  -  S)/(r  -  S)]  ^  ^  exp  [-vS  (r  -q^  )]  (13) 

where 


V 


M- 

-2 

sec 


2 

p/KS 

earth's  gravitational  constant 


-6  3 

1.  5362  X  10  (earth  radii) 


COMPARISON  OF  POWER  AND  EXPONENTIAL  LAW  FITS 
TO  ATMOSPHERIC  DENSITY  DATA  FROM  SEVERAL  SOURCES 

If  (7)  is  to  be  fitted  to  a  given  atmospheric  model  at  a  point  _2 
(pQ,  q^)  in  a  height  band  of  interest,  an  integral  exponent  n  =  1  +  fi 

may  be  determined  by  an  approximate  linear  fit  to  the  scale  height  data. 
Values  of  n  determined  in  this  way  (see  Fig.  2)  are  given  in  Table  2  for 
several  atmospheric  models  for  the  height  bands  noted. 

In  principle,  the  intercept  with  the  r-axis  of  the  line  (3)  fitted  to 
the  scale  height  data  is  the  value  of  s  to  associate  with  each  n.  If  s  is 
determined  in  this  way,  however,  it  is  sensitive  to  uncertainties  in  the 
values  of  H,  particularly  for  the  smaller  values  of  fi  .  For  this  reason 
a  different  procedure  was  adopted  here.  The  air  density  function  (7)  was 
first  identified  with  each  model  at  a  reference  level  q^,  near  the  mid¬ 
point  of  the  height  band  for  which  n  had  been  determined.  An  average 
value  of  s  was  computed  from  the  two  values  obtained  by  equating  the 
function  to  the  air  densities  of  the  model  at  the  end  points  of  the  band. 

The  results  are  shown  in  Table  2. 

For  any  given  height  band,  the  quality  of  the  fit  of  an  air  density 
distribution  such  as  (7)  to  an  atmospheric  model  may  be  expressed  in 
terms  of  the  ratios  of  the  respective  air  densitites.  Calculations  of  this 
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TABLE  2.  VALUES  OF  n  AND  s  IN  SELECTED  HEIGHT  BANDS  FOR 
VARIOUS  ATMOSPHERIC  MODELS 


Model 

Height-Band 

n 

s 

(km) 

A 

150-200 

5 

0.  9993 

200-400 

11 

0.  9580 

B 

150-200 

4 

1.  0135 

200-400 

6 

1. 0010 

C 

150-200 

4 

1.  0064 

200-400 

9 

0.  9716 

D* 

150-180 

2 

1.  0203 

220-400 

5 

1. 0026 

*  The  scale  height  in  the  Paetzold  standard  model  is  roughly 

constant  in  the  interval  180-220  km,  and  the  corresponding 

air  density  distribution  may  be  represented  either  by  (7a)  or 

a  power  function  adapted  to  this  interval. 

type,  based  upon  the  data  in  Table  2,  have  been  made  and  the  results  are 
given  in  Tables  3  and  4  for  the  reference  points  180  and  300  km,  respec¬ 
tively.  It  can  be  seen  that  the  power  laws,  in  general,  represent  their 
respective  models  to  within  a  few  parts  per  hundred  over  the  height  band 
of  the  fit.  Thus,  the  reference  point  (p^,  q^)  may  be  chosen  anywhere 

within  the  height  bands  studied  without  materially  changing  the  quality  of 
the  fit  over  the  band. 

The  fits  of  the  power  laws  to  the  Nicolet  data  and  the  several  pro¬ 
posed  model  atmospheres  may  be  compared  to  the  fits  obtained  with  the 
exponential  law  (11)  by  referring  to  Tables  5  and  6.  These  data  were  cal¬ 
culated  on  the  basis  of  the  published  values  of  scale  height  at  the  reference 
levels  indicated  in  the  tables.  It  is  apparent  that  for  height  bands  of  any 
appreciable  width  the  power  functions  with  integral  exponents  give  closer 
approximations  to  the  models  than  do  the  exponential  functions.  Several 
illustrative  examples  are  shown  in  Figs.  3  and  4. 
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0.  1716 

0. 3212 

0.  4581 

0,  5926 

0.  7008 

0. 9591 

1. 0000 

0.  9013 

*  The  density  scale  heights  used  for  the  calculations  appearing  in  this  table  were  deduced  from  the 
relationship  H  =  (1  +  /3  )  H'  [Nicolet,  1961]. 
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0.  0007 

< 

Q. 

ro 

350  km 
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0.  7362 

0.  9263 

1. 0000 

1.  0364 
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1. 0000 

1. 0471 

0. 9701 

0.  8632 

007 
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0.  7  348 

0.  5803 
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0.  9976 

1. 0000 
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0.  8454 
0.7233 

0. 5865 

0.  3150 

0.  1668 

0.  0853 

Height 
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<km) 
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Fig,  3:  Comparison  of  Exponential  and  Power  Function  Air  Density  Dis¬ 
tributions  Fitted  to  Models  A  and  B . 
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RELATIONSHIP  BETWEEN  a  AND  e  DERIVED  FROM  A  MODIFIED 
POWER  LAW  REPRESENTATION  OF  THE  ATMOSPHERIC  DENSITY 


Several  important  previous  studies  have  been  carried  out  in  which 
the  vertical  distribution  of  air  density  in  the  integrals  (1)  and  (2)  has 
been  represented  by  the  exponential  form  (11)  [e.  g.  -  King-Hele,  Cook 
and  Walker,  I960;  Parkyn,  I960].  In  this  paper,  a  modified  power  func¬ 
tion  with  integral  exponent  will  be  used  to  represent  the  air  density  p. 

In  addition  to  giving  generally  closer  fits  to  the  model  atmospheres  than 
do  the  exponential  distributions,  such  power  functions  also  permit  a  more 
ready  development  of  the  integrals  (1)  and  (2).  The  algebraic  forms 
which  result  from  the  integration  of  (1)  and  (2)  depend  in  a  simple  way 
upon  the  value  of  the  exponent  and  may  easily  be  tabulated  for  various 
choices  of  n  [Peirce,  1929].  These  forms  should  find  application  in  com¬ 
puting  atmospheric  densities  and  comparing  proposed  model  atmospheres 
with  air  densities  deduced  from  observations. 


In  the  treatment  to  follow,  a  simple  series  relationship  (21)  between 
a  and  e  will  be  obtained  from  (1)  and  (2)  using  the  air  density  function  p'^^. 
The  development  will  contain  the  parameters  n  and  s  and  will  be  carried 
out  by  integrating,  with  suitable  approximations,  the  ratio  A  a/  Ae.  The 
later  evaluation  of  (21)  will  be  performed  for  n  =  6,  s  =  1  for  which  p^^' 
will  give  an  acceptable  representation  of  the  shape  of  the  Kallmann  model 
from  150  to  400  km  (See  Table  4).  These  values  have  been  chosen  in 
order  to  present  certain  preliminary  results.  In  general,  the  parameters 
p^,  n  and  s  are  to  be  adjusted  to  fit  the  power  function  to  the  best  avail¬ 
able  atmospheric  density  data. 

The  definite  integral  in  (1)  may  be  written 

(i)  2ir  2  2  3/2  -2  (i) 

I  -  f  (1  -  e  cos  E)  (1  -  e  cos  E)  p  dE, 

0 

i  =  0,  1,  2,  3,  4.  (14) 


For  values  of  eccentricity  greater  than  0.  02,  the  major  contributions  to 
the  integrals  I'^'  occur  at  or  near  perigee  q^.  This  suggests  the  use  of 
the  approximation 


(i)  (i)*  2  3/2 

I  =  (1  -  e  ) 


Zn  -2  (i) 

/  (1  -  e  cos  E)  p  dE. 

0 


(15) 


As  an  illustration  (14)  and  (15)  have  been  programmed  for  the  IBM  7090 
computer  and  evaluated  for  n  =  6,  s  =  1  for  typical  values  of  the  para¬ 
meters  e  and  q.  The  computed  values  of  l(^^  and  I^^'  ,  given  in  Table  7 
in  units  of  p  ,  permit  a  comparison  of  the  sensitivity  of  the  integrals  to 
the  choice  of  air  density  function  and  to  the  approximation  noted. 
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The  substitution  of 


(1)  -n  _ 

p  =  C(n,  s)r(r-s)  ;  C(n,  s)  =  pQ(q^  -  s)  / 


r  =  a(l  -  e  cos  E) 

2  2  1/2  2  1/2 
(1  -  e  cos  E)  »  (1  -  e  ) 

into  (1)  and  (2)  leads  to  the  equations 


3  2  1/2  -n 

Aa  =  -  6  C  a  (1  -  e  )  (a  -  s)  [F(  qin)  +  e  G(Ti;n)] 

2  2  3/2  -n 

Ae  =  -  6  C  a  (1  -  e  )  (a  -  s)  G(  T|;n) 


(16) 

(17) 


where 

■n  =  ae  (a  -  s)  ,  0  =  q  =  1 

ZiT  -n  -(2n  -  1)  /2 

F(q;n)  =  (1  -  q  cos  E)  dE  =  2iT(l-q‘‘)  f{q;n) 

2it  -n  2  "(^^  ■  l)/2 

G(q:n)  =  /q  cos  E(  1  -  q  cos  E)  dE  =  2ir  (1  -  q  )  k(q;n). 

The  functions  f  and  k  are  polynomials  in  q  whose  coefficients  are  readily 
deduced  from  tables  [Peirce,  1929].  The  order  of  f  is  n  -  2,  n  even; 
n  -  1,  n  odd.  The  order  of  f  is  n  -  1,  n  even;  n  -  2,  n  odd.  Explicitly. 

for  2  S  n  g  11, 


n- 1 


n- 1 


5  6 


35 


f(q;n)  -  1  +  2  C^_3q  +  8  16  ^n-?’^  ■'■l28^n-9’’ 


8  _y  n-1 
^  256  n-11 ^ 


nq  1  n-2  2  1  n-2  4  5  n-2 

k(q:n)  ="y  |_1  +  h  +~C_  q  +  rrC, 


•4 ^8^n-6’l  ■'■64^n-8’l  ''■128  ^n-io’l 


7  n-2  8l 

C_  ,„q  J 


where  the  C^_p  ,j=l,  2;p  =  3,  4,  11,  are  the  binomial  coefficients. 

Although  the  use  of  p^^^  or  p^^^  in  (1)  and  (2)  will  lead  to  relationships 
similar  to  (16)  and  (17),  the  particular  forms  of  (16)  and  (17)  are  conven¬ 
ient  for  the  following  development. 

The  quotient  Aa/Ae  may  be  treated  as  a  derivative  to  provide  the 
first  order,  nonlinear  equation 
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+  .f/k)/(l  -A 


(18) 


The  product  e  (f/k)  is  small  compared  to  unity  for  the  entire  lifetime  of 
the  satellite.  This  suggests  that  the  function  f/k  be  expanded  in  a  Taylor 
series  in  powers  of  tj  -  t]q  ,  where  ti^  is  the  initial  value  of  t].  In  turn, 

T)  —  may  be  expanded  in  powers  of  e  -  e^  under  the  approximation 

T]  q  ^/(‘Iq  -  s  +  e).  When  these  expansions  are  introduced  into  (18)  it 
may  be  reduced  to 

1  ^  2  3 

d  In  a  _ 

de 


=  A  (l  +  e)+A,  {e-e)+A  (e-  e^)  +  A  (e  -  e  ) 

0  1  o'  2  o'  3  0 


3  4 

+  e  +  e  +  0  (e  ). 


(19) 


If  the  symbols  f^^^  =  i'f^  and  k^^^  =  ilk^,  i  =  0,  1,  2,  3,  designate  the 
functions  f  and  k  and  their  first  three  derivatives  evaluated  at  the 
constants  Aq,  A  ,  A-,  and  A.  may  be  written 


> 

o 

II 

o 

(a) 

*I  = ’'/'"o  •  ‘ V 

(b) 

(20) 

*2  =  ^2  •  ''l>''‘>0  ■  “  * 

(c) 

3 

A3  =  (y^  -  Zy^  +  y3)/(q^  '  "  +  %> 

(d) 

where 


y  = 

0  0  0 


(q,.  -  8) 


^1  “  q  -  s  +  e 
0  0 

"o  '“o"’ 

- - 2 

<8o'‘*V 


^1  ^0  '  ^0 


f„(k  -kk)-fkk+fk 
0  '  1  0  Z'  1  0  2  2  0 


qo  (q^'s)  fo(2kQkjk^-k^  -kp  k3)+fj(kpkj  -k^j  ^/Vo 

5  0 
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The  exponential  function  obtained  from  the  integration  of  (19)  may 
be  expanded.  The  resultant  series  may  then  be  truncated  and  collected 
in  the  form 

4  . 

a  =  (e  -  e^)^  (21) 


j  =  0 


where  the  Bj  are  functions  of  the  initial  values  a^,  ^through  the  relations 


o 

l( 

(a) 

“  *0  'o  *0  *0  ■ 

(b) 

^2 

=(ao/2)  (1  t  a/  +  Aj  t  4 

'0  -  ^  '0'' 

(c) 

^3 

=(ao/6)  (Aq^  +  3AqA^  +  5A^ 

f  2  A^  -  2  A,  4„) 

(d) 

=(aQ/12)  (4  Aq  A^  +  3  A3)  . 

(e) 

(22) 


In  the  preceding  development  it  has  been  tacitly  assumed  that  the 
earth's  atmosphere  is  spherically  symmetrical.  For  the  purpose  of  de¬ 
veloping  a  satellite  orbit  theory  in  the  presence  of  air  drag  for  perigee 
heights  below  300  km  a  spheroidally  symmetrical  model  of  the  earth's 
atmosphere  is  to  be  preferred.  If  the  spherically  symmetric  atmospheric 
model  is  referenced  to  the  earth's  equator,  a  latitude  dependent  factor 
exp  (  -  €  sin^  §/  H')  [King-Hele,  Cook  and  Walker,  I960]  can  be  intro¬ 

duced  to  take  account  <?f  the  earth's  oblateness.  Here  §  is  the  geocentric 
latitude  and  c  the  earth's  ellipticity  (<  csi  0.003353).  Modified  in  this  way 
the  air  density  function  may  be  written 

*)  =  p(qQ,  0)  (r/q^)  [(q^  -  s)/(r  -  s)]%xp  (-q^e  sin  §/H|j). 

2 

The  factor  exp  ("q^  «  sin  $  Hq)  may  be  expanded  in  terms  of  the  small 

parameter  €  /H^  <0.6  and  the  theory  developed  in  terms  of  a  solution 
due  to  a  spherical  atmosphere  plus  perturbative  terms  arising  from  the 
oblateness . 

The  angle  §  is  related  to  the  orbital  inclination  i  (essentially  con¬ 
stant),  the  argument  of  perigee  w  and  the  true  anomaly 

V  =  cos  ^  [(cos  E  -  e)/(l  -  e  cos  E)] 

2  2  2 

through  the  relation  sin  $  =  sin  i  sin  (v  +  u).  Thus,  the  oblateness 
effect  depends  on  the  parameters  c,  q^,  i,  u,  e  and  E,  although  in 

(1)  and  (2),  the  integrations  are  carried  out  under  the  assumption  that 
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only  E  varies  during  a  single  revolution  from  perigee  to  perigee.  Con¬ 
sequently  £he  gross  influence  of  oblateness  may  be  conveniently  inter¬ 
preted  in  terms  of  "external"  contributions,  independent  of  E,  due  to  the 
precession  of  the  perigeal  line  over  an  oblate  earth  and  "internal"  con¬ 
tributions,  dependent  on  E,  due  to  a  single  orbital  passage  of  the  satellite 
over  an  oblate  earth.  For  orbits  of  moderate  eccentricity  (e>  0.  02)  both 
the  "external"  and  "internal"  contributions  are  suppressed  in  ratios  such 
as  (18).  The  "external"  contributions  cancel  directly  and  the  "internal" 
contributions  effectively  cancel  because  most  of  the  drag  occurs  at  or 
very  near  perigee.  Hence  expression  (21)  is  relatively  insensitive  to 
oblateness  effects  for  eccentricities  e>0.02  [King-Hele,  Cook  and 
Walker,  I960]. 


DISCUSSION  OF  RESULTS 

Equation  (21)  has  been  programmed  for  the  IBM  7090  computer  and 
calculations 'of  a  (n  =  6;  s  =  1),  based  upon  orbital  data  compiled  by  Walker 
[Walker,  1961],  have  been  made  for  all  of  the  satellites  of  interest 
launched  through  I960.  Certain  representative  results  are  shown  in 
Tables  8  and  9.  Values  of  a  computed  from  (21)  and  the  corresponding 
values  computed  from  the  K-H  C  W  [King-Hele,  Cook  and  Walker]  theory 
for  an  oblate  atmosphere  are  tabulated.  Also  given  are  values  of  a  de¬ 
duced  from  observations  of  the  periods  of  the  satellites  [Walker,  1961]. 
Values  of  a  computed  from  (21)  using  Smithsonian  mean  eccentricities  for 
1958  Y  (Izsak,  I960)  are  given  in  Table  9  along  with  the  counterpart 
Smithsonian  data. 

The  agreement  of  the  predictions  of  (21)  with  the  K-H  C  W  theory, 
with  the  Smithsonian  mean  elements  and  with  the  observed  values  is  quite 
good  over  the  range  of  eccentricities  (e>  0.  02)  for  which  (21)  is  expected 
to  hold.  The  differences  in  the  lengths  of  the  semimajor  axis  given  either 
by  the  K-H  C  W  theory  or  by  Smithsonian  data  and  those  given  by  (21)  are 
less  than  5  km  with  a  single  exception.  The  exception,  a  difference  of  13 
km,  occurs  for  1957/3  for  the  small  value  of  eccentricity  e  =  0.  022  at 
which  (21)  might  be  expected  to  yield  poorer  results,  particularly  for 
orbits  with  a  near-critical  inclination  angle  (see  below).  General  agree¬ 
ment  of  a  similar  quality  between  the  values  computed  from  (21)  and  the 
observed  values  may  also  be  noted. 

No  simple  well  ordered  analysis  of  the  effects  of  atmospheric  oblate¬ 
ness  is  possible  since  they  depend  onamedleyofparameters.  There  should 
be,  however,  a  strong  dependence  on  both  the  orbital  inclination  and  the 
argument  of  perigee  due  to  the  appearance  of  sin^tj)  in  the  latitude -depend¬ 
ent  factor  [King-Hele,  Cook  and  Walker,  I960].  The  close  agreement  of 
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COMPARISON  OF  VALUES  OF  a  COMPUTED  FROM  (21)  WITH  PREVIOUS  THEORY  AND 
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the  data  from  (21)  with  the  K-H  C  W  theory  for  satellites  such  as  1957  fi  , 
1958  6  and  1959  which  have  significantly  different  values  of  the  oblate¬ 
ness  parameters  illustrates  the  fact  that  atmospheric  oblateness  has 
little  influence  on  the  variation  of  a  with  e  for  e>0,  02.  The  influence  of 
atmospheric  oblateness  becomes  more  important  for  the  less  eccentric 
orbits  particularly  for  orbits  vith  inclination  near  the  critical  angle 
(63.4*)  where  the  precessional  rate  of  the  argument  of  perigee  is  small, 
as  with  1957/3  and  1958  6^.  If  u  sweeps  out  one  or  more  revolutions,  as 
with  1958  €,  an  averaging  effect  is  present  which  tends  to  make  the  oblate¬ 
ness  effect  small.  This  is  illustrated  by  the  fact  that  the  differences  in  a 
between  the  values  given  by  (21)  and  the  K-H  C  W  theory  is  greater  for  the 
satellites  1957>/3  and  1958  with  e  =  0.022  and  e  =  0.017,  respectively, 
than  for  the  satellite  1958  «  with  e  =0.011. 
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